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Background Short-range correlations (SRC) are an important ingredient of the dynamics of nuclei. 

Purpose An approximate method to quantify the magnitude of the two-nucleon (2N) and three-nucleon (3N) short-range 
correlations (SRC) and their mass dependence is proposed. 

Method The proposed method relies on the concept of the "universality" or "local nuclear character" of the SRC. We quantify 
the SRC by computing the number of independent-particle model (IPM) nucleon pairs and triples which reveal beyond- 
mean-field behavior. It is argued that those can be identified by counting the number of nucleon pairs and triples in a 
zero relative orbital momentum state. A method to determine the quantum numbers of pairs and triples in an arbitrary 
mean-field basis is outlined. 

Results The mass dependence of the 2N and 3N SRC is studied. The predictions are compared to measurements. This 
includes the ratio of the inclusive inelastic electron scattering cross sections of nuclei to 2 H and 3 He at large values of 
the Bjorken variable. Corrections stemming from the center-of-mass motion of the pairs are estimated. 

Conclusions We find that the relative probability per nucleon for 2N and 3N SRC has a soft dependence with mass number 
A and that the proton-neutron 2N SRC outnumber the proton-proton (neutron-neutron) 2N SRC. A linear relationship 
between the magnitude of the EMC effect and the predicted number of proton-neutron SRC pairs is observed. This 
provides support for the role of local nuclear dynamics on the EMC effect. 

PACS numbers: 25.30.Fj,24.10.-i,13.60.Hb 



I. INTRODUCTION 

We define the nuclear packing factor (NPF) as the 
fraction of the nuclear volume that is occupied by nu- 
cleons. A rough order of magnitude estimate of the 
NPF can be arrived at using uniform spheres for the 
nuclear and nucleon density. The nuclear radius Ra 
can be reasonably determined from Ra = 1.2(fm) A 1 / 3 . 
It is not obvious what value of the nucleon radius 
should be used. In models of relativistic heavy-ion col- 
lisions it is customary p] to use expulsion distances d, 
which simulate the hard-core NN repulsion, of the or- 
der of 1 fm, corresponding with rjv ~0.5 fm. This leads 
to NPF=0.07. A recent reanalysis of electron scattering 
data resulted in a root-mean-square charge radius of the 

proton r c p = y / ( r 2)=0.897(18) fm 2\. Assuming that 

the r° v is an estimate of the proton and neutron radius 

/ r"(fm) \ 3 

one arrives at NPF= ( J =0.42. It is clear that the 

computed NPF is very sensitive to the adopted value of 
the nucleon radius. The estimate of the NPF on the ba- 
sis of Tp should be considered as an upper limit. Indeed, 
the established value of the nuclear saturation density of 
0.17 nucleons/fm corresponds with a mean internucleon 
distance of 1.8 fm implying that rjy < 0.9 fm. 

From the above, it is clear that one expects that the 
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nucleus is more like a saturated quantum liquid than a 
gas of freely moving nucleons. Accordingly, the nuclear 
wave functions receive large corrections from short-range 
(SRC) and long-range correlations. These days it is com- 
mon practice to implement the effect of SRC in nuclear 
computations. Examples include the calculations of ma- 
trix elements for double-/3 decay [3], of event simulations 
in heavy- ion collisions [J], and of hadron transparencies 
in nuclei [5]. 

The EMC effect [5] is the reduction of the cross section 
for leptonic scattering off a nucleon bound in a nucleus 
relative to that of a free nucleon (mass Mm)- The EMC 
effect was observed in Deep Inelastic Scattering (DIS) 
experiments on nuclei at high virtual-photon virtualities 
Q 2 = q 2 - uj 2 > 2 GeV 2 for Bjorken x B = ^m^Z in 
the range 0.3 < xb < 0.7. The ratio of per nucleon 
cross sections is denoted by R — 4tt where a A is the 
cross section for leptonic scattering from the target A. 
The magnitude of the EMC effect can be quantified by 
means of the slope — ^p- [7]. Another remarkable fea- 
ture of the ratio R is that it adopts a constant value 
(this factor is commonly referred to as the SRC scaling 
factor ci2(A/D)) for 1.5 < ig < 2 and moderate values of 
Q 2 [SHIP]. It has been suggested [TI] that the a 2 (A/D) 
can be related to the high-momentum components of the 
nuclear wave functions. A phenomenological linear re- 
lationship between the a^i^AjU) and the magnitude of 
the EMC effect expressed as — Jjp- has been observed 
[121114) . This indicates that the magnitude of the Euro- 
pean Muon Collaboration (EMC) effect may be driven 
by SRC. In this picture the magnitude of the EMC effect 
is (partly) related to the temporal local density fluctu- 
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ations which are induced by the high virtualities of the 
leptonic probe. Recent measurements 7 5 corroborate this 
relation between the local nuclear environment and the 
magnitude of the EMC effect. 

Given an arbitrary nucleus A(N, Z) we address the is- 
sue of quantifying the number of two-nucleon (2N) pairs 
prone to SRC and the number of 3N triples prone to 
SRC. Along the same lines we investigate to what extent 
the mass dependence of the NN SRC can be captured by 
some approximate principles. We wish to develop a ro- 
bust method which is applicable to any nucleus from He 
to Pb. From this method we expect, for example, that 
it allows one to study the mass dependence of the SRC 
without combining results from various types of calcula- 
tions. 

Momentum distributions contain the information 
about IN, 2N, 3N, . . . properties of the nuclear ground 
state. Over the years various methods to compute 
the nuclear IN and 2N momentum distributions have 
been developed. Ab-initio calculations which solve the 
Schrodinger equation with realistic nucleon-nucleon in- 
teractions are available for light nuclei like 4 He [T51 - 
IT7] . For medium- weight nuclei (12 < A < 40) trunca- 
tion schemes based on cluster expansions can be adopted 
[15] , Correlated-basis function theory has been applied 
to compute ground-state densities and momentum distri- 
butions for doubly-closed-shell nuclei from 12 C to 208 Pb 
[T51 [213] , Thanks to the enormous progress in theoret- 
ical many-body nuclear physics and the availability of 
nuclear momentum distributions in a broad mass range, 
times are ripe to learn more about SRC, for example by 
mapping its A and isospin dependence. It remains no- 
toriously difficult, though, to establish quantitative rela- 
tionships between observables and the computed momen- 
tum distributions [HI I2D - I25] . Here, we do not attempt 
a high-precision calculation of momentum distributions. 
Our goal is to gather insight into the mass and isospin 
dependence of the SRC from stylized facts of momentum 
distributions. 

In a mean-field model fluctuations are completely ig- 
nored. The SRC induce spatio-temporal fluctuations 
from the mean-field predictions for the nuclear density 
distributions for example. As a result of SRC, realistic 
nuclear wave functions reflect the coexistence of single 
nucleon (mean-field) structures and cluster structures. 
The clusters account for beyond mean-field behavior. As 
the nucleon-nucleon interaction is short ranged, the clus- 
ters attributed to SRC are predominantly 2N. The cen- 
tral result of this paper asserts that the amount of 2N 
and 3N SRC in nuclei can be reasonably quantified by 
counting the number of nucleon pairs and triples in a zero 
relative orbital state in a mean-field ground-state wave 
function. In order to quantify the isospin dependence of 
the 2N and 3N correlations, additional information about 
the spin dependence of the clusters is necessary. 

This paper is organized as follows. Sect.[lTJis devoted 
to a discussion of momentum distributions and of how 
they can be used to quantify the mass and isospin depen- 



dence of SRC. In Sect. IIIII we address the issue whether 
inclusive electron scattering data can be linked to the 
number of correlated 2N and 3N clusters. Thereby, we 
deal with both the a2(A/D) coefficient and the magni- 
tude of the EMC effect. 



II. QUANTIFYING NUCLEAR 
CORRELATIONS 



In this section we start from stylized facts of nuclear 
momentum distributions in order to arrive at criteria to 
quantify the 2N and 3N SRC in nuclei. Our focus is on 
their mass dependence. 



A. Nuclear momentum distributions 

In this subsection we provide the definitions and nor- 
malization conventions of the nuclear momentum distri- 
butions used here. For the sake of the simplicity of the 
notations, we will only consider the positional degrees-of- 
freedom. Unless stated otherwise the spin - and isospin 
degrees-of-freedom are not explicitly written in the ex- 
pressions. 

The one-body momentum distribution of nuclei is de- 
fined as 



Pi 



(*) =7^3 / d ^ J dr(^-(^) Pl {T U ri) , (1) 



where pi(ri,r/) is the one-body non-diagonal density 
matrix 



Pi(ri,f{) = J {df2-N}^A(ri,r2,f 3 ,...,f A ) 

x^> A (f;,r 2 ,f 3 ,...,r A ) . (2) 

Here, ^>a is the ground-state wave function of the nucleus 
A and the notation 

{dri- N } = dndr l+1 . . . df A , (3) 

has been introduced. For (*f? A \ ^a) = 1, one has that 



dkP l (k = 1 



(4) 



We introduce relative and center-of-mass (cm.) co- 
ordinates of nucleon pairs in coordinate ^r*i2,i?i2^ and 

momentum space I fci2 , P12 



ri2 



r\ - r 2 
V2 



fci - k 2 

~7T 



Rr 2 - (5) 



r> kl + fc2 ra\ 
P 12 = — , (6) 



V2 
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and define the two-body momentum distribution in the 
standard fashion as 



(fcl2,A 2 ) 



f 



dr 12 / dR 12 / dr{ 2 / dRl 



xe «fel2 ' (fl2-r{ 2 ) e iP 12 ■ (R 12 -R{ 2 ) 

X P2 {?i2,Ri2]r( 2 ,Ri^j . (7) 



Here, p 2 \?X2, R12', r{ 2l R{ 2 ) is the two-body non- 
diagonal density matrix 



P2 



[ri 2 ,Ri 2 \f{ 2 ,R{ 2 ^j = p 2 



R{. 



^ +ri 2 + Ri 2 _ -fi2+Ru +r{ 2 

{dr 3 -N} *a (n,r 2 ,f 3 , . . . , r A ) $ A (r{, f 2 , r 3 , . . . , r A ) 

I 



2 zfl 



-,r 2 = 



^12 + #12 



72 



(8) 



One has the normalization condition 



dk 12 / dP 12 P 2 (k 12 ,P 12 ) = 1 



(9) 

In a spherically symmetric system, the two-body momen- 
tum distribution P 2 (k\ 2 , P\ 2 ^j depends on three indepen- 
dent variables. One of the most obvious choices [5^ is 



(1 



fcl2 1,1 Pl2 



&12-Pl2 



(10) 



where ^% 12 p 12 is the angle between P\ 2 and k\ 2 . 

The distributions P\ (kj and P 2 \ki2,Pi2j reflect all 
information about one-nucleon and two-nucleon proper- 
ties contained in the ground-state wave function. Other 
quantities can be directly related to them. Here, we list 
some of the most frequently used ones. 

The two-body cm. momentum distribution is defined 

as (dP l2 = P? 2 dP 12 dn Pl2 



W12) = / dk 12 / dn Pl2 P 2 (k 12 ,P 12 



(11) 



The quantity P\ 2 P 2 (P12) dP\ 2 is related to the probabil- 
ity of finding a nucleon pair in A with cm. momentum 
P\2 =| P\ 2 I irrespective of the value and direction of the 
relative momentum k\ 2 of the pair. The P 2 (P12) receives 
contributions from the proton-proton, neutron-neutron, 
and proton-neutron pairs 



P2 (Pi. 



ppp 

r 2 



(Pr. 



r>nn 
r 2 



(P12) 



pr (P12) 



(12) 



In a spherically symmetric nucleus, it is convenient to 
introduce the quantities 

ni(fe) = JdClkPt^tj , (13) 

n 2 (k 12 ,P 12 ) = J dfl kl2 J dQ Pl2 P 2 (£12, A2) -(14) 



The quantity n\ (k) k 2 dk gives the probability of finding 
a nucleon with a momentum in the interval [fc, k + dk\. 
The n 2 {k\ 2l P\ 2 )k\ 2 dk\ 2 Pl 2 dP\ 2 is the combined prob- 
ability of finding a nucleon pair with a relative mo- 
mentum in [k±2, ki 2 + dki 2 ] and cm. momentum in 
[P 12 ,P l2 + dP l2 ]. 



B. Mean-field approximation and beyond 

A time-honored method to account for the effect of 
correlations in classical and quantum systems is the in- 
troduction of correlation functions. Realistic nuclear 
wave functions | , 5) can be computed after applying a 
many-body correlation operator to a Slater determinant 

I tyMF\ 



1 



A 



Q\Q I $f 



g 1 *r > . (15) 



The nuclear correlation operator Q is complicated but as 
far as the short-range correlations are concerned, it is 
dominated by the central, tensor and spin correlations 



S 



II ( 1 ~ 9c(nj) + ftr(rij)SijTi ■ fj 



(16) 



where g c (ri 2 ), ftr(ri2), fariru) are the central, tensor, 
and spin-isospin correlation function, S the symmetriza- 
tion operator and S\ 2 the tensor operator 



Si 2 = -o-^i • r 12 <x 2 ■ r 12 - a l ■ 3 2 



xy 2Ml (n n2 ) Pi ® s 2 \ 2 _ M 



24tt 



M L 



(17) 
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The operator S12 admixes relative two-nucleon states 
of different orbital angular momentum, is operative on 
triplet spin states only, and conserves the total angular 
momentum of the pair. 

We stress that the correlation functions cannot be con- 
sidered as universal and that in some many-body ap- 
proaches, particularly for light nuclei, they do not ap- 
pear. The momentum distributions which result from 
the calculations depend on the interplay between many 
factors, including the choices made with regard to the 
nucleon-nucleon interaction, the single-particle basis (if 
applicable), the many-body approximation scheme, .... 
As a matter of fact, different nucleon-nucleon interactions 
and many-body approaches may produce, particularly in 
the region of SRC (short distances/high momenta), mo- 
mentum distributions which are very similar (see, e.g. 

Refs. nana esi es). 

The g c {f 12) quantifies how strongly two point-like nu- 
cleons treated as quasi-particles, are spatially correlated 
when they are a distance r\2 apart. The g c (ri2) gives rise 
to local density fluctuations about the mean-field predic- 
tions from the reference state | ). The GD g c {ri2) 
(computed for nuclear matter) from Fig. [I] is not very 
different from the one for monoatomic molecules in a liq- 
uid. Indeed, for r 12 — > one has that the GD g c (r 12 ) — > 1 
which reflects the fact that nucleons have a finite size (or, 
in other words they are subject to a nucleon-nucleon in- 
teraction with a hard core). For values of r\2 which are 
larger than a few times the diameter of a nucleon, the 
gc ( r i 2 ) 0. From this we conclude that the fluctua- 
tions from the MF densities are confined to short inter- 
nucleon distances. Therefore, the 2N SRC are a highly 
local property and are insensitive to the properties of 
the other surrounding nucleons. This is the fundamental 
reason why SRC can be considered as "universal" [T7] , 
Whereas a large model dependence for the g c is observed, 
the f tT seems to be much better constrained. We have 
added the squared D-wave component of the deuteron 
wave function ^fjj (fci 2 ) in Fig. Ill Obviously, the mo- 
mentum dependence of | ft T (^12) | 2 and the deuteron 
momentum distribution nu =| ^ u (&12) | 2 are highly 
similar. 



Eqs. and Q can be written as 

Pi (k) = P<°> (k) + pW (k) , (18) 

ft(*ia,^ia) = P 2 m (*12,A 2 ) +P 2 {1) (*ia,Aa) ■ 

(19) 

The P^ and P% are the mean-field parts and are fully 
determined by the Slater determinant | X S I ^ F ). After 
inserting the expressions ([!]) and |2| into the Eq. ( 13 ) 
one obtains 



4 0) (fc) = J <in k pw (fe) = 1 J2 (?h + 1) s, 



n h lhih 



nhlhjh 



X drr 3i h {kr)il) nh i hjh (r) 



(20) 



where ji(r) is the spherical Bessel function of the first 
kind and the sum extends over all occupied single-particle 
states. The < S nh i h j h < 1 is the occupation probability 
of the corresponding single-particle state. The presence 
of short-range and long-range correlations leads to occu- 
pation probabilities smaller than one. With the adopted 
normalization convention of Eq. Q one typically obtains 
that 



dkk 2 n { ° ] (k) 



0.6-0.8 



(21) 



or, about 60—80% of the nucleons are mean-field like. We 
stress that a considerable fraction of this depletion can 
be attributed to long-range correlations, an effect which 
is not considered here. 

The distribution k 2 n^\k) as it can be computed 



from Eq. ( 20 1 is reminiscent for a phenomenon which 
is confined to a certain scale, or, in other words, it is 
Gaussian like. The typical scale is determined by the 
Fermi momentum kp ps 250 MeV. This is illustrated 
in Fig. [2] where we show the momentum dependence 

of the k 2 n { °\k) for 12 C, 56 Fe and 208 Pb as computed 
with Woods-Saxon (WS) wave functions. For the sake of 

curiosity we have fitted the computed k 2 n^\k) with a 
Boltzmann distribution 



The effect of the correlation functions on the momen- 
tum distributions can be roughly estimated from their 
squared Fourier transforms. The effect of the tensor 
correlation function is largest for moderate relative mo- 
menta (100 < ki2 < 500) MeV. For very large ku, the 
g c is the dominant contribution. The harder the <7 c (ri2) 
the stronger the effect of correlations. We stress that 
in the plane-wave impulse approximation, the SRC con- 
tribution to the (e, e'pp) cross section is proportional to 
| 5c (fci 2 )| 2 E5]. 



After introducing the wave functions of Eq. (151, 



the one-body and two-body momentum distributions of 



4-7T 



(2TrM N kT) 



3/2 



k 2 exp - 



A- 2 



2M N kT 



(22) 



The results of the one-parameter fit are shown in Fig. [2] 
The fit is remarkably good for Carbon and gets increas- 
ingly inaccurate with increasing mass number. From the 
fit of the Boltzmann distribution we obtain kT ps 12 MeV 
(C), kT fa 14 MeV (Fe), kT ps 16 MeV (Pb). Accord- 
ingly, for the IPM part of the momentum distribution, 
the typical energy exchange per momentum degree-of- 
freedom \kT is of the order of 6-8 MeV. 

The correlated part k 2 n^(k\ on the other hand, is 
reminiscent of the nucleus as a system of interdependent 
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FIG. 1: (color online). The radial and momentum dependence of a central and some tensor correlation functions. The central 
correlation function "GD" is for nuclear matter and from Ref. 1 30 . The tensor correlation function "Pieper" is for 16 and 
from Ref. [28], the "CBF" one is for 16 and from Ref. [2j5|, and the "cluster" one is for ie O and from Ref. [31], * D (k 12 ) is 
the I = 2 component of the non-relativistic deuteron wave function generated with the Paris potential [321 133| (not to scale). 



nucleons and is obviously non-Gaussian. In contrast to 
the mean-field part nf^ , the correlated part extends 
over "all" momentum scales. Or, in other words the 2N, 
3N, . . . correlations generate a fat momentum tail to the 
n\ (fc). The high momentum tails to n\ (fc) have a very 
similar form for all nuclei, including the deuteron, which 
alludes to some universal character of SRC [T7] . 

It has been theoretically predicted [3"3r[3"rj] and exper- 
imentally confirmed in semi-exclusive A(e, e'p) measure- 
ments [37] that the major fraction of the (k > hp) 
strength is contained in very specific parts of the single- 
nucleon removal energy-momentum phase space, namely 
those where the ejected nucleon is part of a pair with 
high relative and small cm. momentum. This is the so- 
called ridge in the spectral function [35] which reflects 
the fact that high-momentum nucleons in the one-body 
momentum distribution are related to 2N dynamics with 
two nucleons which are close and move back-to-back with 
approximately equal and opposite momenta. 

From recent calculations |26| of the two-body momen- 



tum distributions in 3 He and 4 He the following conclu- 
sions could be drawn. At high relative momenta and 
small cm. momenta, the cm. and relative motion of the 
pair is decoupled, an effect which is reminiscent of 2N 
SRC. For the correlated pn pairs the relative motion can 
be described by the high-momentum part of the deuteron 
wave function. This suggests the following expression for 
the correlated part of the pn two-body momentum dis- 
tribution 

n { 2 ] (2&f < fc i2 , P 12 < 150 MeV) 

^a pn (A,Z)n D (k 12 )FP n (P 12 ) , (23) 

where a pn (A, Z) is a proportionality factor related to 
the number of correlated proton-neutron pairs in the nu- 
cleus A Z relative to the deuteron and nn (^12) is the 
high-/ci2 deuterium momentum distribution. Further, 
the F pn (P12) is the cm. distribution of the correlated 
pn pairs. It corresponds with that part of P 2 (P12) 
of Eq. that stems from pn pairs with a zero rel- 
ative orbital angular momentum li 2 = and a total 
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FIG. 2: (color online). The computed k 2 n { °\k) versus k for 
the nuclei 12 C, 56 Fe and 208 Pb and a Boltzmann fit including 
the error bars. The extracted values of kT are 12.0±0.5 MeV 
(C), 14±1 MeV (Fe), and I6±l MeV (Pb). The calculations 
are performed with WS single-particle states. The adopted 
normalization convention is J dk k 2 nf\k) — 1. 



spin S = 1. The proposed scaling behavior (23) can 
be attributed to the dominance of the tensor correla- 
tions at medium relative momenta and the fact that 
I ftr (&12 > kp) | v E'_D(fci2)| , two qualitative obser- 
vations which can made from Fig. [I] 



C. Quantifying two-nucleon correlations 

We suggest that the significance of 2N correlations in 
a nucleus A(N, Z) is proportional to the number of rel- 
ative ?i2 = states [21]. There are experimental re- 
sults supporting this conjecture. First, in high-resolution 
16 0(e, e'pp) 14 N measurements performed at the electron 
accelerators in Amsterdam [55] and Mainz [35] , the quan- 
tum numbers of the target nucleus and the residual nu- 
cleus are unambiguously determined. For the transitions 
to low-lying states in the residual nucleus, the eightfold 
differential cross section for the exclusive (e, e'pp) reac- 
tion has been studied as a function of the initial cm. 
momentum P 12 of the proton-proton pair which is in- 
volved in the reaction process. This has provided insight 
into the quantum numbers of the pairs involved in the 



12 



(ri 2 ),A 12 (£«)) 



reaction process. We denote by 

the orbital wave function corresponding with the relative 
and cm. motion of a nucleon pair. For the ground-state 
(g.s.) to g.s. transition, for example, 



16 



O(0+,g.s.) + e 



C{0 + ,g.s.)+e'+pp, (24) 



the active diproton resides in a state with quantum num- 
bers \li2 = 0, A 12 = 0) at lower P 12 and |Z l2 = 1, A 12 = 1) 
at higher P 12 . Two independent calculations from 
the Pavia and Ghent groups have demonstrated that 
the largest contributions from SRC to the eight-fold 
cross section are confined to low P 12 values [39 . This 
provides direct evidence of pp correlations being con- 
fined to \l±2 = 0, Ai 2 = 0) pairs. In that sense, the 
16 0(e, e'pp) 14 N measurements nicely confirmed the back- 
to-back picture of SRC: diprotons are subject to SRC 
whenever they happen to be close (or, in a relative Zi 2 = 
state) and moving back-back (or, in a state with P 12 « 
which corresponds with A 12 = 0). 

High-resolution (e, e'pn) measurements which have the 
potential to access the pn correlations are very challeng- 
ing [10]. Theoretical (e,e'pn) calculations [5rJl |4"T1 |4"2"] 
have predicted that the tensor parts of the SRC are re- 
sponsible for the fact that the correlated pn strength is 
typically a factor of 10 bigger than the correlated pp 
strength. Calculations indicated that the tensor corre- 
lations are strongest for pn pairs pairs with "deuteron- 
like" \li2 = 0,S = 1) relative states gHH3- Recently, the 
dominance of the pn correlations over pp and nn ones has 
been experimentally confirmed [4"31 144) . 

Accordingly, a reasonable estimate of the amount of 
correlated nucleon pairs in A(N, Z) is provided by the 
number of pairs in a lyi = state. In order to deter- 
mine that number for a given set of single-particle states, 
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one needs a coordinate transformation from (ri,r 2 ) to 



T12 



V2 



',Ri 



ri+r 2 \ 
V2 J' 



For a harmonic oscillator 



(HO) Hamiltonian this transformation can be done with 
the aid of Moshinsky brackets [15] 

I mil (fi)naia {f 2 )]LM L ) = ^ 

(nuh2NuAu;L \ ni*in 2 / 2 ;£) 

X I "12*12 (^12) ^Vi 2 A 12 (ti 12 y,LM L ). (25) 

We define the interchange operator for the spatial, 
spin, and isospin coordinate as 



N 
N 



I 0.1 









pp(T=l) 


s=o — ■— 




n 12 =0,S=0 x 





10 

mass number A 



100 



V12 = V12 (ri,f 2 )Vi2 (cti, a 2 )Vi2 (n,f 2 ) 



(26) 



After introducing the spin and isospin degrees-of- 
frccdom, in a HO basis a normalized and antisym- 
metrized two-nucleon state reads (o^ = (miiji^i)) 

' 1 2 ' ' > ™° V 2 (1 + ^1^2) 
x (1-Pi 2 )|ai (fi)aa (r 2 );JM) 

- E E E E E 



LM^ n 12 ii2 N12A12 SAf s TA/ T V 2 0- + ^"102) 

x [1 - 

x(ni2/i 2 A r i 2 Ai2;i | nihn 2 l 2 ;L) 

f*i h L) 

xjJ 2 LS ; ; S (LM L SM S | JM) 
[h 32 J) 

y.{ l -h l -t 2 \TM T ) 



"42*12 (?"l2),A^i 2 Ai2 (-R12) 



LM L ,SM S ,TM T 



(27) 



where we have used the shorthand notation j = \J2j + 1. 

With the above conventions one has that the total 
amount of proton-neutron pairs can be obtained from 
a sum over all pn pairs in the nuclear ground state 



E E E na 

JM ai < a !, d!<aj 



aia 2 ; JM \otia 2 ; JM) na = NZ , 

(28) 

where a p F and ctp denote the Fermi level for the proton 
and neutron. Similar expressions hold for the number of 
proton-proton and neutron-neutron pairs 

^r^E E E 

JM ai <a p F a 2 <a p F 

na {aia 2 ; JM \ ai a 2 ; JM) na , (29) 

JM ai<a™ «2 

na (a 2 a 2 ; JM |aia 2 ; JM)„ a . (30) 
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FIG. 3: (color online). The computed values for Z ^_ 1 ~ ) N PP , 

jv(jv-i) N nn , and Npn(S) which represent the predicted 
fraction of the pairs which are prone to SRC. The re- 
sults are obtained for HO single-particle wave functions with 
hbj(MeV) = 45. A" i - 25. A~ § and for the target nuclei 4 He, 
9 Be, 12 C, le O, 27 A1, 40 Ca, 48 Ca, 56 Fe, 63 Cu, 108 Ag, and 
197 Au. 



Starting from the Eq. (27) one can compute in a HO 
single-particle basis how much a pair wave function with 
quantum numbers 

"12*12 (ri 2 ),iV 12 A 12 (i? 12 )] LM L ,SM S ,TM T ) (31) 



contributes to the sum- rules of Eqs. ( |28|29|30[ ). This can 
also be done for any other basis \nljm) of non-relativistic 



single-particle states. In that case, the adopted proce- 
dure involves an extra expansion of \nljm) in a HO basis 

\nljm) = {lmi-m a \ jm)ip n ij(r)Y lmi (n)x| m . 

mi Tn a 

= E (/ *n*K B i( r Wrai{r)\ \n H ljm) , (32) 

where <pn H i{T) are the radial HO wave functions. A two- 
nucleon state can then be expressed in a HO basis for 



which the Eq. (27 1 can be used to determine the weight 



of the pair wave functions of Eq. (31) 



The IPM pp pairs are mainly subject to the central 
SRC which requires them to be close. This implies that a 
reasonable estimate of the number of IPM pp pairs which 
receive substantial corrections from the SRC is given by 
an expression of the type 



N PP (A,Z)=J2 E E 



na (aia 2 ; JM\ V^~° \a x a 2 ] JM) 7 



(33) 



where V^ 2 is a projection operator for two-nucleon rel- 



ative states with li 2 — 0. A similar expression to Eq. ( 33 ) 



holds for the nn pairs. For the pn pairs it is important to 
discriminate between the triplet and singlet spin states 



N pn{s) (A,Z)=J2 E E 

JM ai<a p F a2<aj- 

na ( ai a 2 ; JM\ V l £=°V§ Ka 2 ; JM) r 



(34) 



In Fig. [3] we display some computed results for the 
Npp, N nn , and N pn ^ for 11 nuclei. The selection of 
the nuclei is motivated by the availability of inclusive 
electron-scattering data and covers the full mass range 
from Helium to Gold. We have opted to display the 
results relative to the sum rule values of the Eqs. (28) 
and ( 29 ) , which allows one to interpret the results in 



terms of probabilities: given an arbitrary pair wave func- 
tion, what is the chance that it has zero orbital rela- 
tive momentum and a specific spin quantum number. 
In a naive IPM picture for 4 He, the pp pair is in a 
\lx2 = 0, S = 0,T = 1) state. As this 2N configuration is 
prone to central SRC effects, the corresponding probabil- 
ity is 1. The physical interpretation is that for 4 He "all" 
IPM pp-pair wave function combinations receive correc- 
tions from SRC. For a medium-heavy nucleus like 56 Fe or 



(>:; 



Cu we find 



A, 



0.1, which leads one to conclude 



that about 90% of the IPM pp pair wave functions do 
not receive corrections from central SRC. For the heav- 



iest nucleus considered here (Au) 



2N„ 



= 0.06, which 



Z(Z-l) 

means that only about 190 out of the 3081 possible pp 
pair combinations are subject to SRC. 

Comparing the mass dependence of the pp and nn re- 
sults of Fig. [3] one observes similar trends. For the pn 
results a softer decrease with increasing A is predicted. 
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/ 

/ ri2(ni2/l2) 
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FIG. 4: (color online). Use of Jacobi coordinates in the ppn 
system. 



There are about three times as many pn(T = 0) states 
than pn(T = 1) states with I12 = 0. This would be trivial 
in a system with only spin and isospin degrees of freedom. 
In a system in which the kinetic energy plays a role and 
in which there are spin-orbit couplings, we cannot see 
any trivial reason why this should be the case. In this 
respect, we wish to stress that for most nuclei discussed 
N Z. A stronger criterion for selecting nucleon pairs at 
close proximity is imposing n\ 2 = in addition to l± 2 = 
and we have added also those results to Fig. [3j We find 
the results of Fig. [3] robust in that the A dependence and 
magnitudes are not very sensitive to the choices made 
with regard to the single-particle wave functions. All the 
results of Fig. [3] are displayed on a log-log plot and can 
be reasonably fitted with a straight line, pointing towards 
a power-law mass dependence A a for the N pp , N pn and 



MS)- 



D. Quantifying three-nucleon correlations 

In order to quantify the magnitude of the 3N corre- 
lations for an arbitrary A(N, Z), we build on a recent 
paper by Feldmeier et al. [I7j. There, it is pointed out 
that 3N correlations can be induced without introducing 
genuine three-body forces. In terms of the correlation 
operators of Sect. |IIB] three-body correlations will natu- 
rally emerge in cluster expansions by means of operators 
of the type g c (r 12 ) g c (na), g c (n.2) ftr (ru), 

The strongest source of three-nucleon correlations is 
the tensor correlation operator acting on the (S — 1,T = 
0) channel of the pn states of Eq. (27). We consider 
ppn configurations and explain one possible mechanism 
to create a correlated state (see also Fig. [4| . In the un- 
correlated wave function one has a n(l)p(2) pair in a 
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\li2 — S12 = 1 ; Ti2 = 0) state and a p(2)p(3) pair in a 
\hz = S 23 = ; T 2 3 = 1) state. Accordingly, both pairs 
are in relative I — states. In Ref. [T7] it is explained that 
these two pairs can be brought into a correlated three- 
nucleon status by flipping the spin of proton 2. In the 
correlated part of the wave function one has an n(l)p(2) 
pair in a | Z 12 = 2 Si 2 = 1 ; Tyx = 0) and an p(2)p(3) pair 
in a I Z23 = 1 5*23 = 1 ;T 2 3 = 1) state. This configura- 
tion can be energetically favorable through the presence 
of the strong tensor correlation in the pn pair. Indeed, 
the energy gain through the tensor induced n(l)p(2) cor- 
relation can compensate for the energy loss of breaking 
the pairing in the p(2)p(3) pair. 

Given A(N, Z) we propose to find all the antisym- 
metrized 3N states with orbital quantum numbers 



(ni 



0Zl2 



0) 



(35) 



"(12)3 = ^(12)3 

in the IPM wave function and identify them as the dom- 
inant contributors to 3N SRCs. This corresponds with 
seeking for those 3N wave-function components where 
all three nucleons are "close". This can be techni- 
cally achieved by constructing antisymmetrized 3N states 
starting from a MF Slater determinant, and performing 
a transformation from the particle coordinates (fx, f 2 , f 3 ) 

to the internal Jacobi coordinates m.2, ^(12)3, ^123 



R12 ~ V2f 3 
73 ' 



Rx23 — 



V2R 12 



'■3 



(36) 



One readily finds for uncoupled three-nucleon states in a 
HO basis [45] 

\n 1 l 1 m h (fx) , n 2 l 2 m h (r 2 ) , n 3 l 3 m h (f 3 )) = 

E E E E E 

LM L n 12 li2 N 12 A 12 L 1 M Ll n(i 2 )3'(i2)3 

x E E E 

N123A123 m ll2 m Al2 ™i (12)3 ™A 123 

x(lxmij 2 m l2 I LM L )(l 12 m h2 A 12 M Al2 \ LM L ) 
x(A 12 M Al2 l c mi c I LxM Ll ) 
x{l {12)3 m lll2)a A 123 M Al23 I LiM Ll ) 
x(ni 2 h 2 N 12 A 12 ;L \ nihn 2 h;L) 

x ("(12)3'(12)3^123Ai23;il | Ni 2 Ai 2 n 3 l 3 ; Li)p 

X \n 12 l 12 m ll2 (f 12 )) I ri(12)3^(12)3^ ( i 2) 3 ( f (12)3)) 

X |^Vl23Al23M Al23 (iW)) , (37) 

where we have adopted the notation (. . . | . . .}p for the 
Standard Transformation Brackets (STB) [45] . 

Antisymmetrized (a) uncoupled three-nucleon states 
can be obtained from the three-nucleon wave functions 



of Eq. ( 37 1 using the interchange operator of Eq. ( 26 ) 



\a a m a ,a b m b ,a c m c ) a = 

[1 \a a m a (fi),a b m b (f 2 ),a c m c (f 3 )) 

[1 - V\ 2 \ \a b m b (fx) , a c m c (f 2 ) , a a m a (f 3 )) 

[1 - V12] \a c m c (fx) , a a m a (f 2 ) , a b m b (f 3 )) .(38) 




FIG. 5: (color online). The mass dependence of 
the amount of ppn triples with quantum numbers 
\ni2 = 1 12 = 0,71(12)3 = Z(i2)3 = 0). The results can be 
reasonably fitted with a power law O^SA 1 ' 58 * 20 . The re- 
sults are obtained for HO single-particle wave functions with 
fko(MeV) = 45.A~i - 25.A"§ and for the nuclei 4 He, 9 Be, 
12 C, le O, 40 Ca, 48 Ca, 56 Fe, 63 Cu, 108 Ag, and 197 Au. 



The total number of ppn triples can now be expressed as 
Z(Z-l) 



N- 



E E 

a a ,a b <a F „ a a <a™ m a m b m c 



E 



na (a a m a ,a b m b ,a c m c \a a m a a b m b a c m c ) na ,(39) 

which allows for a stringent test of the analytical deriva- 
tions and their numerical implementation. Along similar 
lines to those used to derive the number of correlated 2N 



clusters in Eq. ( 33 1 , the number of ppn triples with the 



orbital quantum numbers of Eq. ( 35 ) can be obtained 
from 



N ppn (A,Z) 



E E E 

a a ,a b <a p „ a c <a™ m a m b m c 



I I T)ni2=0,Zl2=0T-,™(12)3— 0,2(12)3— 



r (12)3 



\a a m a , a b m b , a c m c ) 



(40) 



We associate the N ppn (A,Z) with the number of ppn 
SRC triples. The A dependence of N ppn (A, Z) is dis- 
played in Fig. [5] There is striking linear correlation be- 
tween the logarithm of the mass number and the loga- 
rithm of the number ppn triples which are close in the 
MF ground-state wave function. 



III. RESULTS 

In this section we discuss how our predictions for the 
number of correlated 2N pairs and correlated 3N triples 
can be connected with experimental results from inclu- 
sive electron scattering. 



10 



Separation of the correlation and mean-field 
contributions 



We start with illustrating that the separation of the 
mean-field and correlated contributions to the inclusive 
A(e, e') cross sections is feasible. In order to achieve this, 
we use stylized features of the n\ (k) in Monte-Carlo (MC) 
simulations to illustrate that a separation between the 
mean-field n^\k) and the correlated n^\k) part can be 
made in the A(e,e') signal. 

We assume that quasi-elastic single-nucleon knockout 
e + A — > e' + (A — 1) + N is the major source of A(e, e') 

strength. With q(u),q), pa(M a ,0), pa-i(Ea-i,PA-i), 
Pf(E]y,pf) we denote the four-momenta of the virtual 
photon, of the target nucleus, of the residual A — 1 sys- 
tem, and of the ejected nucleon. From energy-momentum 
conservation 



q + PA -PA-i =Pf, 



(41) 



one can deduce for A = 2 a relation between the min- 
imum of the missing momentum p m — pf — q and the 
Bjorken scaling variable xb for fixed Q 2 [35]. The results 
are shown in Fig. [6] Obviously, for Q 2 > 1.5 GeV 2 and 
xb > 1.5 one mainly probes nucleons with a momentum 
well above the Fermi momentum for the deuteron. For 
finite nuclei the situation is more involving as A — 1 rep- 
resents an additional degree of freedom which can carry a 
fraction of the transferred four-momentum. We have per- 
formed MC simulations for a fixed energy of the imping- 
ing electron beam e, and a fixed electron scattering angle 
6 e . The p m for a mean- field nucleon is drawn from the 
MF part nf\k) of n 1 (k). For a correlated nucleon the 



p m is drawn from n x (fc). Parameterizations for n^''{k) 

,0) 



(0), 



and n\ (k) are obtained from [35] 



4 0) 
.0) 



(fc) = A^e 
(k) = A (1) e 



■[l + 0(k 2 )}, 



(42) 



C^e~ D{1)k \ (43) 



where A^>, B^°\ A^\ CO) and £>0) depend on A. 

In Fig. [7j we compare the xb distribution of simula- 
tions for the mean-field and correlated pa rt of one nu- 
cleon knockout in 12 C. As stated in Eq. (21 ), the number 
of events is normalized as J dk k 2 rif\k) — 0.7. For 
xb > 1.5, the events originate almost uniquely from 



,0) 



(k). 



B. Two-body correlations 



Following the experimental observation ]§\ [TOl S3 that 
the ratio of the inclusive electron scattering cross sections 
from a target nucleus A and from the deuteron D 



a A (x B ,Q 2 
ct d {x B ,Q 2 



(44) 
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FIG. 6: (Color online). Relation between the minimum of the 
missing momentum for the deuteron and is at various 

values of the frour-momentum transfer Q 2 . 
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FIG. 7: (Color online). Simulation of one-nucleon knockout in 
12 C with separated mean field and correlated momentum dis- 



tribution. The number of events is normalized as in Eq. (21 1 
The electron kinematics is determined by e. 
9e = 18°. 



5.766 GeV and 



scales for 1.5 < Xb ^$ 2 and moderate Q 2 , it has been 
suggested [37] to parameterize the a A in the following 
form 



° A (1-5 < x B < 2, Q 2 ) = -as (A/D) a 2 (A, x B ,Q 2 ) , 

(45) 

where a 2 (A, xb,Q 2 ) is the effective cross section for scat- 
tering from a correlated 2N cluster in nucleus A. Assum- 
ing that CT2 is some local function which does not depend 
on the target nucleus A 



a 2 (A,x B ,Q 2 ) wct 2 (A = 2,x B ,Q 2 ) ^ a D (x B ,Q 2 



(46) 
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one can rewrite Eq. ( 45 ) as 



In this simplified reaction-model picture, which ignores 
amongst other things the cm. motion of pairs in finite 
nuclei, the quantity — a 2 (A/D) can be connected with 
the number of correlated pairs in the nucleus A. As- 
suming that all pn pairs contribute one would expect 
that for the relative amount of correlated two-nucleon 
clusters a 2 (A/D) ~ A. Based on the observed domi- 
nance of correlated pn pairs over pp and nn pairs |43| . 
and the universality of the deuteron-like high momentum 
tail of the correlated two-body momentum distribution 
(23), we suggest that the correlated pn pairs contribut- 
ing to the a%{A/D), are predominantly (T = 0, S = 1) 



correlated pn(5 = 1) pairs dominate the signal [13] 2H 
[52]. 

As pointed out in Refs. [TTJ [25], the cross section for 
the exclusive (e, e' NN) reaction can be written in a fac- 
torized form as 

a A (e, e'NN) = KF NN (P 12 )a eNN (k 12 ) , (48) 

where P\ 2 {ki 2 ) is the cm. (relative) momentum of the 
correlated pair on which the photoabsorption takes place 
and if is a kinematic factor. The above expression is 
valid in the plane-wave and spectator approximation for 
electron scattering on a pair with zero relative orbital mo- 
mentum. The (J e NN stands for the elementary cross sec- 
tion for electron scattering from a correlated 2N pair with 
relative momentum fc 12 . The <J e jvAr contains the Fourier- 
transformed correlation functions g c (ki 2 ) and ft T (k\2). 



pairs and that 02 (A/D) is proportional to the quantity An analytic expression for a epp is contained in Ref. 



N, 



(A,Z) defined in Eq. (34 1 



pra(S=l) 

In Ref. [3H] the ratio of Eq. (|44|) has been calculated 



with spectral functions obtained from state-of-the-art nu- 
clear matter calculations in the local density approxima- 
tion for the correlated part and A(e, e'p) scattering data 
for the mean-field part [1§1 [50]. The calculations sug- 
gested large FSI effects, whilst the plane-wave calcula- 
tions did not exhibit the scaling present in the data at 
1-5 < Xb- In Refs. [TTJ [55] it is argued that a complete 
treatment of FSI in this kinematics needs to include in- 
elastic channels in the rescattering and that this cancels 
part of the elastic FSI contribution included in Ref. [48] . 
The results in Ref. [35] seem to refute the validity of 



Eq. (45), which hinges on the assumption that the FSI 



effects on correlated pairs in a nucleus are almost identi- 
cal to those in the deuteron in a high-momentum state. 
This requires that for 1.5 < the FSI is primarily in the 
correlated pair and that the remaining A — 2 nucleons act 
as spectators. Such a behavior is suggested by the calcu- 
lation of the quasi-elastic cross sections in Ref. [51] and 
by a space-time analysis of the nuclear FSI at xb > 1 
carried out in Ref. |11) where it is stressed that the rein- 
teraction distances are < 1 fm, supporting the idea that 
the first rescattering should be very similar to FSI in the 
deuteron (see a recent discussion in Ref. [H]). Therefore 
the assumption of Eq. ( 45 1 seems a reasonable one for 



light nuclei where the amount of rescatterings is of the 
order of 1. For medium-heavy and heavy nuclei, the aver- 
age amount of rescatterings is larger than 1 and it has to 
be verified if the assumption still holds. The settlement 
and clarification of all the cited issues related to the role 
of FSI in inclusive reactions requires further studies with 
a full reaction model. 

In a finite nucleus correlated pairs can have a non- 
zero cm. momentum. This cm. motion is a correction 
factor when connecting the measured a2(A/D) to the 
number of correlated pn pairs N pn /s—i) (A, Z) . We aim to 
provide an estimate for this correction factor. Therefore, 
we consider the two-nucleon knockout reaction e + A — > 
e' + (A— 2)+N+N following the break-up of a correlated 
2N cluster. For an inclusive cross section, the tensor 



and has been tested against data in Ref. [53]. 

As argued above, in order to link the exclusive cross 
section of Eq. ( 48 1 to the inclusive ones contained in the 
Eq. ( 45 1 one assumes that a epn s» a e jj and one introduces 
a proportionality factor ^Vp n (s=i)(A, Z) which counts the 
number of correlated pn pairs in A. With the scaling 



relation of Eq. (48) for the (e,e'pri) reaction, one can 



transform the ratio of Eq. ( 47 1 into a form which accounts 



for the cm. motion of the pair 



a 2 (A/D) = 



A 



J PS dk 1 2dP 12 N pnis=1) (A,Z)FP n (P 12 )a eD (fc 12 ) 



■jN pn(s=1) 



J PS dki 2 a eD {ki 2 ) 
(A,Z) [ dP 12 F pn (P 12 ) , 

JPS 



(49) 



where the integrations extend over those part of the cm. 
momentum phase (PS) included in the data. A basic 
assumption underlying the above equation is that the 
factorization of Eq. (23) approximately holds. The com- 



puted widths of the cm. momentum distributions for the 
correlated pn pairs contained in Table 1 indicate that the 
major fraction of the pn pairs has P\ 2 < 150 MeV which 
is within the ranges for the validity of Eq. ( 23 ) . 



In line with our assumption that the correlated pairs 
are dominated by pn in a relative 3 Si state, F pn (P12) 
can be expressed as the conditional two-body cm. mo- 
mentum distribution 



F pn (Pi 2 )=P% n (PufSt) 
Figure M shows calculations for the P% " (P) 

pP™ I plTc. ' 



(50) 
and 

P2' 1 (P| iSi) for 12 C. The cm. distribution of correlated 
pn pairs (P pn (Pi 2 )) can be well parameterized in terms of 
a Gaussian distribution. The widths <r c m obtained from 
a Gaussian fit to P|" 1 (Pi2| 3 <Si) are given in Table [ij 

To estimate the cm. correction factor we have per- 
formed MC simulations of pn knockout with and with- 
out inclusion of the cm. motion. This amounts to 
drawing the cm. momentum from F pn (Pi 2 ) = 5(P\ 2 ) 
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FIG. 8: (Color online) The pn cm. distribution Pi 2 Pj n (P12) 
and PfaPf 1 (Pi 2 | 3 Si) = P? 2 F pn (P 12 ) for 12 C. The calcula- 
tions are performed with HO single-particle states and adopt 
the normalization convention J dP\i P^P^™ (P12) = NZ. 



and F pn (P 12 ) ~ e 2 "c.m. , where <r c . m . is the A depen- 
dent width. For 1.5 < xb < 2, the initial momen- 
tum distribution of the correlated pair is given by cor- 
related part of the tw o-b ody momentum distribution 
\ki2,Pi2). The Eq. (23 1 states that the n 2 1] (k 12l P 12 ) 



,(!)/ 



can be considered universal. As illustrated in Fig. [T] one 
has Tin (^12) ~ |/tr(^i2)| 2 - As the relative momentum 
distribution is approximately proportional to the tensor 
correlation function, we draw k\2 from the distribution 
^ ; i2 1 ftr (^12) 1 2 • Energy conservation reads 



(q + PA- Pa-2 - Pa) 2 = p) 



(51) 



where q(ui,q), pa(M a ,0) and pa-2(Ea-2, ~(p s + Pm)) 
are the four-momenta of the virtual photon, target nu- 
cleus and residual A — 2 system, respectively. The vir- 
tual photon interacts with one of the nucleons, result- 
ing in a fast nucleon pf(Ef,pf) with Pf = p m + q and 
a slow nucleon p s (E s ,p s ). With the aid of Eq. (51 ), one 
can calculate the a^-distribution of the simulated events. 
We apply the kinematics of the Jefferson Lab ( JLab) ex- 
periment E02-019 [TO]: e t = 5.766 GeV and 9 e = 18°. 
The average < Q 2 > of the generated events (including 
cm. motion) in the xb region of interest is 2.7 GeV 2 . 
This value, which is A-independent, agrees with the one 
quoted in Ref. [10] . 

The results of our simulations are summarized in 
Figs. [9] and [10] Fig. [9] shows the xb — kvx scatter plot 
of 10 simulated events with and without inclusion of 
cm. motion for 12 C. In both situations the mass differ- 
ence between inital and final state causes a small shift to 
lower xb compared to the deuteron case. Second, we ob- 
serve considerable shifts in the distribution of the events 



in the (k\2,XB) plane due to cm. motion. In Fig. 10 
one can observe how cm. motion considerably increases 
the number of events with 1.5 < xb < 2. The impact 





A 


(Jem. 


cm. correction factor 




L2 C 
56 Fe 

208p b 


115 MeV 
128 MeV 
141 MeV 


1.64 ±0.23 

1.70 ±0.27 

1.71 ±0.29 



TABLE I: The second column gives the width of the cm. 
distribution of correlated pn pairs. The third column provides 
the computed cm. correction factor. The errors represent the 
dependence on the choice of correlation function. 



of the cm. corrections increases with growing xb- Ex- 
perimentally, the a 2 (A/D) coefficient is determined by 
integrating data for 1.5 < xb < 1.85. We estimate the 
cm. correction factor by the ratio 

# simulated events with inclusion of cm. motion 
# simulated events without inclusion of cm. motion 

in this xb region. The resulting correction factor for 
several nuclei is contained in Table [I] We performed the 
simulations with the three different correlation functions 
f tT in Fig. [T] The dependence of the result on the choice 
of correlation function is represented by the error of the 
cm. correction factor. 

Fig. [3] quantifies the fraction of all possible pn pairs 
which are prone to SRC relative to the total amount of 
possible pn pair combinations. In our picture one has 
N pn (s=i) = 1 f° r D. This means that we do interpret 
the I12 — component of the deuteron wave function as 
the IPM part which receives large corrections from tensor 
SRC. The per nucleon probability for a pn SRC relative 
to the deuterium can be defined as 



N pn{s=1) (A, Z) 



N + Z N, 



pra(S=l) 



(A = 2,Z = 1) 



.4 



N pn(s=1) (A, Z) 



Similar expressions hold for the per nucleon pp SRC and 
the per nucleon nn SRC 



;Npp(s=o)(A, Z) 



N 



N nn(S=Q){A, Z) . 



(54) 



The re sult s of the per nucleon probabilities are collected 
2 H, 



in Fig. 11 Relative to H, the per nucleon probability 
of pn SRC are 2.20, 3.63, 4.73 times larger for Carbon, 
Iron, Gold. Along similar lines, relative to the "free" pp 
system the per nucleon probability of pp SRC are 1.39, 
2.34, 3.11 times larger for Carbon, Iron, Gold. 

In Fig. [12] we compare our predictions computed with 
the aid of the Eq. (49) with the extracted values of 



a,2(A/D). We have opted to correct the predicted a 2 co- 
efficients and not the data for cm. motion. We stress 
that the cm. correction factor cannot be computed in a 
model-independent fashion. For light nuclei our predic- 
tions tend to underestimate the measured a 2 . This could 
be attributed to the lack of long-range clustering effects 
in the adopted wave functions. Indeed, it was pointed 
out in Ref. [54] that the high-density cluster components 
in the wave functions are an important source of correla- 
tion effects beyond the mean-field approach. For heavy 
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FIG. 10: (Color online) Histogram of the xb distribution of 



12 C(e,e'pn) MC simulations with {F(P 12 ) ~ e 2<7 c. m . ) and 
without (F(Pi2) ~ <5(Pi2)) inclusion of cm. motion. The 
kinematics is the one of the JLab experiment E02-019 [10] : 
<a = 5.766 GeV and 8 e = 18°. 



nuclei our predictions for the relative SRC probability 
per nucleon do not saturate as much as the data seem 
to indicate. In Ref. [TT] the authors estimated the mass 
dependence of ai by means of an expression of the type 
a2 ~ J d 3 rp^ lF (r). Using Skyrme Hartree-Fock densities 
Pmf(^) a power-law of A 012 emerged. After normaliz- 
ing to the measured value for a 2 ( 12 C/D) the predicted 
power-low dependence agrees nicely with the data. 

We stress that final-state interactions (FSI) represent 
another source of corrections which may induce an ad- 
ditional A-dependent correction to the data. FSI of the 
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FIG. 11: The mass dependence of the per nucleon probability 
for pn SRC relative to the deuterium. 



outgoing nucleons with the residual spectator nucleons, 
could shift part of the signal strength out of the cuts ap- 
plied to the experimental phase space (or likewise move 
strength in) and decrease (or increase) the measured 
cross section and the corresponding a-i coefficient. 

In Fig.[l3]we display the magnitude of the EMC effect, 
quantified by means of ~4^- versus our predictions for 



the quantity ^N pn is—i-\ or, the "per nucleon probability 
for pn SRC relative to the deuteron" . We stress that the 
numbers which one finds on the x-axis are the results of 
parameter- free calculations. We consider the "per nu- 
cleon probability for pn SRC relative to the deuteron" as 
a measure for the magnitude of the proton-neutron SRC 
in a given nucleus. Obviously, there is a nice linear rela- 
tionship between the quantity which we propose as a per 
nucleon measure for the magnitude of the SRC and the 
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FIG. 12: (Color online) The computed values for the a^A/D) 
for various nuclei. The data are from Refs. [9] 1101 |14j . The 
shaded region is the prediction after correcting the computed 
values of a,2(A/D) for the cm. motion of the pair. The cor- 
rection factor are determined by linear interpolation of the 
factors listed in Table U The width of the shaded area is de- 
termined by the error of the cm. correction factors. 



Sec. IIIB| it has been suggested [3] to parameterize the 
inclusive A(e,e!) cross section in the following form 



~aa (A/ 3 He) 



03 (xb,Q 



<j A (2.25 < x B < 3,Q 2 ) 

where 03 (xb, Q 2 ) is the cross section for scattering from 
a correlated 3N cluster which is once again assumed to 
be A independent. Inserting Eq. (56 1 into Eq. (55), one 
obtains 



(57) 

Notice that in the kinematic regime where 3N correla- 
tions are expected to dominate (2.25 < Xb) the experi- 
mental situation is unsettled. For example, the recently 
measured a3( 4 He/ 3 He) ratios [TO] are significantly larger 
than those reported in Ref. [5]. 

Similar to the per nucleon pn SRC of Eq. (53) we 



define the per nucleon probability for a ppn SRC relative 
to 3 He as 





0.6 




0.55 




0.5 


QD 


0.45 


X 
T3 


0.4 


cc 


0.35 


1 


0.3 




0.25 




0.2 




0.15 




3 4 

(2/A) N pn(S=1) 



FIG. 13: The magnitude of the EMC effect versus the com- 
puted per nucleon number of correlated pn pairs. The data 
are from the analysis presented in Refs. [7] 1141 I55| . The 
fitted line obeys the equation —43. = (0.108 ± 0.028) + 



|iV (s= ir (0-074 ±0.010). 



magnitude of the EMC effect. 



C. Three-body correlations 

The measurements of Refs. [§] [TO] indicate that the 
ratio of the inclusive cross sections 



,r 3 Hc 



(XB,Q 2 ) ' 



(55) 



approximately scales for 2.25 < xb < 3.0. Along sim- 
ilar lines as those used in quantifying the 2N SRC in 



3 N ppn (A,Z) 3 

A N ppn (A = 3,Z = 2) A Nppn {A Z) ' (58) 



where we used the fact that N ppn (A = 3, Z = 2) = 1 in 
our framework. The results of the per nucleon probability 



of ppn SRC are collected in Figure 14 



The quantity of Eq. (58) can be linked to a,3(A/ 3 He) 
under the condition that corrections stemming from cm. 
motion of the correlated ppn triples, FSI effects, . . . are 
small. Under those idealized conditions one would have 



a 3 (A, 3 He)« jN ppn (A,Z) 



(59) 



In the naive assumption that all 3N pairs contribute to 
the a,3(A/ 3 He) ratio, one expects an A 2 dependency. We 
suggest that only ppn triples in a "close" configuration 
contribute and we count the number of SRC triples with 
the aid of the Eq. ( 40 ) . The ppn contributions will be 
larger than the pnn ones due to the magnitude of the 
electromagnetic coupling. Correlated triples should have 
at least one pn pair due to the dominant character of 
the tensor component. In Fig. [14] we show the predic- 
tions for the a3(A/ 3 He) coefficient as computed with the 
Eq. ( 59 1 and compare it to the data. We stress that the 



experimental situation is largely unsettled and that nei- 
ther the data nor the theoretical calculations have been 
corrected for cm. motion and FSI effects. For Helium 
and Carbon our predictions are in line with the experi- 
mental value. For Iron the prediction is about a factor 
of two larger than the experimentally determined ratio 
of cross sections. Our parameter-free calculations repro- 
duce the fact that the mass dependence is much softer 
than the A 2 dependence that one would expect on naive 
grounds. 
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FIG. 14: The mass dependence of the per nucleon probability 
for ppn SRC relative to 3 He. We stress that neither the data 
nor the theoretical calculations have been corrected for cm. 
motion and FSI effects. The data are from Ref. [5]. 



IV. CONCLUSION 

We have provided arguments that the mass depen- 
dence of the magnitude of the NN and NNN correlations 
can be captured by some approximate principles. Our 
method is based on the assumption that correlation op- 
erators generate the correlated part of the nuclear wave 
function from that part of the mean-field wave function 
where two nucleons are "sufficiently close". This trans- 
lates to computing those parts of the two-nucleon and 
three-nucleon wave functions with zero relative orbital 
momentum in order to identify short-range correlated 
pairs and triples. 

We have calculated the number of pn, pp and nn 
I12 = SRC pairs and studied their mass and isospin 
dependence. The A dependence of the magnitude of the 
pp, nn, and pn SRC manifests itself in a power-law de- 
pendence. We found a significant higher per nucleon SRC 
probability for pn pairs than for pp and for nn. To con- 



nect the computed number of SRC pairs to the measured 
<Z2 (j4/D) corrections are in order. Published experimen- 
tal data include the radiation and Coulomb corrections. 
The correction factor stemming from final-state interac- 
tions and from the cm. motion of the correlated pair, 
however, is far from established. We proposed a method 
to estimate the cm. correction factor based on general 
properties of nucleon momentum distributions. Using 
Monte Carlo simulation we find a correction factor of 
about 1.7 ±0.3. Our model calculations for are of the 
right order of magnitude and capture the A-dependcncc 
qualitatively. For small A our predictions underestimate 
the data, while we do not find the same degree of satu- 
ration for high A that the (scarce) data seem to suggest. 

To compute the number of 3N SRC in a nucleus, we 
count the ppn states with three nucleons which are close. 
We have quantified the number of 3N SRC and pro- 
vided predictions for the measured 03 (A/ 3 He) coeffi- 
cients. Our model calculations for the 03 are of the same 
order of magnitude as the (scarce) data but overestimate 
the 56 Fe data point. In this comparison no corrections for 
cm. motion and FSI effects have been made and it re- 
mains to be studied in how far they can blur the connec- 
tion between inclusive electron-scattering data and the 
SRC information. We find a linear relationship between 
the magnitude of the EMC effect and the computed per 
nucleon number of SRC pn pairs. This may indicate that 
the EMC effect is (partly) driven by local nuclear dy- 
namics (fluctuations in the nuclear densities), and that 
the number of pn SRC pairs serves as a measure for the 
magnitude of this effect. 
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